Biomimetic isotropic nanostructures for structural coloration 
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Many species of birds have feathers that are 
brilliantly-colored without the use of pigments. In 
these cases, light of specific wavelengths is selectively 
scattered from nanostructures with variations in in- 
dex of refraction on length scales of the order of vis- 
ible light. ^ This phenomenon is called structural 
color. The most striking examples of structural color 
in nature are iridescent colors created by scattering 
from periodic structures. EJEJ The colors produced 
by these structures change dramatically depending 
on the angle of observation. Nature also produces 
structural colors that have very little angle depen- 
dence. These colors are the result of scattering from 
isotropic structures. SI El El 

In recent years, periodic biological structures have 
provided inspiration for groups trying to make pho- 
tonic materials. Much of this work has 
been motivated by producing a photonic band gap. 
EES] EH EJUS] EH However, Nature's alternative de- 
sign, based on isotropic structures is just starting 
to be explored. Hallam et al have used 
biomimetic random structures to make ultra-thin 
mineral coatings that are brilliant white. Takeoka et 
al recently showed that a wide range of colors with 
very little angle dependence can be produced by mi- 



crogel dispersions. In this communication, we de- 
scribe the self-assembly of biomimetic isotropic films 
which display structural color that is amenable to po- 
tential applications in coatings, cosmetics, and tex- 
tiles. We find that isotropic structures with a char- 
acteristic length-scale comparable to the wavelength 
of visible light can produce structural color when 
wavelength-independent scattering is suppressed. 

We make two types of films that are structurally- 
colored by exploiting the self-assembly of colloidal 
polymer nanoparticles. The first type of sample is a 
thin film on a glass coverslip produced by spin casting 
an aqueous suspension of spheres (Figure |TJl) . While 
monodisperse dispersions form anistropic polycrys- 
talline films, as shown in the inset of Figure a 
mixture of two sizes of spheres ensures an isotropic 
structure as shown in Figure [TJd.EI The blue-green 
film in Figure [I] is made from a bidisperse suspension 
of polystyrene (PS) spheres with mean diameters of 
226 and 271 nm and polydispersity 2% in equal vol- 
ume fractions. 

We image the structure of the film using scanning 
electron microscopy (SEM). The side-view SEM in 
Figure [TJd is of a film comparable to the one in Fig- 
ure [IJi and shows that there is no long range order, it 
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Figure 1: Effect of disorder and order on optical 
properties, a) Photograph of a film of 226 and 271 nm 
PS spheres spin-cast onto an 18 x 18 mm glass cov- 
erslip. Inset: Photograph of a dried sessile droplet. 
b) Side-view SEM image a similar film. The field of 
view is 8.8 jam wide, c) The top set of curves show the 
transmission spectra for the isotropic film pictured in 
(a). The bottom set of curves show the results of the 
same measurement performed with a crystalline sam- 
ple, made by spin-casting 226 nm PS spheres. The 
data were taken with the sample normal to the opti- 
cal path (solid line), at an angle of 30° (dashed red 
line), and at an angle of 80° (dotted blue line), rep- 
resented schematically in the upper left hand corner. 
Inset: Side-view SEM image of a crystalline sample, 
the field of view is 5.3 jam wide. 



Figure 2: Structure of isotropic films, a) SAXS pat- 
tern from an isotropic film, the field of view is 0.15 
nm -1 and the gray values are logarithmic in inten- 
sity, b) The azimuthal average of the experimental 
(red dots) and numerical (blue dashed line) scatter- 
ing pattern. The black line is optical reflectivity data 
taken at an angle of 10° (experiment geometry shown 
in Figure [4| converted to q-space using an effective 
refractive index of n e = 1.24. 

also reveals a representative thickness of 2.3±0.2 /am. 

The optical properties of isotropic and crystalline 
films are quite different, even for samples prepared 
with very similar particles. We compare the trans- 
mission spectra of isotropic and crystalline samples 
in Figure The isotropic film is the same one pic- 
tured in Figure [T^i, the crystalline film is composed 
solely of the 226 nm spheres. Both samples show 
a dip near 500 nm at normal incidence. However, 
when the angle of incidence is changed, the spectral 
position of the dip for the crystalline sample shifts 
and eventually disappears. In contrast, the position 
of the dip for the isotropic sample does not move. 
This illustrates the trade-off in optical performance 
between crystalline and isotropic structures: while 
the crystalline film has more pronounced features at 
some angles, the isotropic film performs consistently 
over a wide range of angles. 

To quantify the structure of the films, we perform 
small-angle X-ray scattering (SAXS) measurements 
on samples composed of spheres with mean diame- 
ters of 226 and 265 nm with 2% polydispersity, spin- 
cast on kapton tape. A typical scattering pattern is 
shown in Figure [2^i. The pattern is dominated by 
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a ring of uniform intensity, indicating that there is a 
well-defined length scale with no preferred direction - 
the structure is isotropic with short-range order. The 
azimuthal average of this pattern is shown in Fig- 
ure [2}d, along with the expected scattering intensity 
from simulations of jammed packings of bidisperse 
spheres with size and number ratios that match our 
experimental system. The position of the first peak 
in I(q) occurs at qo = 0.03 nm" 1 for both simula- 
tion and experiment. The full width at half maxi- 
mum (FWHM) of the first peak, Aq, characterizes 
the range of spatial order £ = 2ir/Aq = 870 nm. In 
powder crystallography, £ describes the crystal do- 
main size, here, £ is only a few particle diameters. 

We directly compare optical reflectivity measure- 
ments to SAXS measurements by plotting the re- 
flected intensity as a function of wavevector q = 
(47rn e /A) cos (0 m /2), where n e is the effective refrac- 
tive index of the material; m is the angle between 
illumination and detection, taking into account re- 
fraction at the film surface; and A is the wavelength 
of light in vacuum (Figure J2Jd). The peaks from both 
measurements match when n e = 1.24. We apply the 
Maxwell-Garnett equation, 
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to calculate the volume fraction of spheres, <p = 
0.46±0.04, where nps is the index of refraction of 
the spheres, taken to be 1.58, and n a i r is taken to be 
1.00. 

Isotropic films can produce structural color with 
little angle dependence, but the film thickness criti- 
cally affects its color, as seen in the inset of Figure [lji. 
Here, we cast a thick film by drying a sessile droplet of 
the same suspension used to make the thin films. In 
thick regions near the center, the film appears white. 
In thin sections near the edge, it appears blue-green. 
This thickness dependence can be understood in the 
following way. The film preferentially scatters wave- 
lengths corresponding to the peak in I(q). In a thin 
film, only these wavelengths will be scattered to the 
detector resulting in a structural color. In a thick 
film, all wavelengths are scattered multiple times and 
reach the detector. 
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Figure 3: Optimizing color by adding absorption, a) 
Photograph of five drop-cast films of 226 and 265 nm 
PS spheres containing carbon black. Sample num- 
bers and [CB] in wt% appear above and below the 
samples, respectively, b) Normalized optical scatter- 
ing spectra recorded at an angle of 20° (experiment 
geometry shown in Figure [4]) for the five samples in 
(a). Inset: Non- normalized scattering spectra for the 
same samples, c) SEM image of the interior region of 
sample 3, the field of view is 7.8 fim wide, a piece of 
CB is indicated with an arrow, d) The width of the 
spectra for each sample in (a) at I(X)/I(\q) = 0.90 



The sensitivity of the color to the film thickness 
requires well-controlled casting procedures which in- 
crease the cost and limit the coated area. There- 
fore, for many potential applications it is necessary to 
eliminate the thickness-dependence. We address this 
by introducing broadband absorption to the bulk of 
the films. The absorption length plays a similar role 
as the thickness: it limits the path length of light 
through the film by absorbing photons that do not 
get scattered within a small distance from the surface. 
Thus, only wavelengths with the strongest scattering 
will escape the film before being absorbed. 
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We make a series of thick films with different ab- 
sorption lengths by drop-casting films with varying 
concentrations of carbon black. The photo in Fig- 
ure [3^i illustrates the effect of adding carbon black 
on the color of the material. From left to right, the 
concentration of carbon black, [CB], is increased. In- 
tuitively, one might expect a mixture of black and 
white to make gray, but here they make blue-green 
over a range of [CB]. The plot in Figure^ shows the 
normalized optical scattering spectra for these five 
samples. As [CB] is increased from 0.02 wt% to 1.80 
wt%, the contrast between scattered intensity at the 
peak and shorter wavelengths is increased, improving 
the color of the sample. When [CB] = 11.2 wt%, the 
scattered intensity is reduced by a factor of ten com- 
pared to the brightest sample (inset in Figure [3]o) , 
and the film appears dark gray. The reflectance peak 
is narrowest for [CB] = 1.80 wt% (Figure |3]i), and 
indeed this sample has a distinct blue-green color. 
Table [I] lists the [CB] for each sample and the cor- 
responding extinction length: the extinction length 
is extrapolated from measurements of aqueous sus- 
pensions of CB. Interestingly, the extinction length 
for sample 3 is 1.3 /im, which is comparable to the 
thickness of the thin film pictured in Figure [IJl 

Table 1: Carbon Black Extinction Length 

Sample [CB] [wt%] Extinction Length [m] 



1 O02 1.3xl0- 4 

2 0.18 1.3xl0- 5 

3 1.80 1.3xl0- 6 

4 3.60 6.5xl0- 7 

5 11.2 2.1xl0~ 7 



The thick films with CB are not iridescent un- 
der omnidirectional illumination, but, under direc- 
tional illumination, the peak wavelength scattered 
does change slightly when the angle between illumi- 
nation and detection is varied. ^ Since the colors 
are the result of single scattering, the position of the 
scattering intensity maximum does not vary with re- 
spect to g, as demonstrated by the spectra in Fig- 
ure H] The red dots connected by a red line in Fig- 
ure represent the azimuthal average of the SAXS 
pattern for a sample with the same [CB] prepared on 
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Figure 4: Color is the result of single scattering. Nor- 
malized reflectance spectra taken over a range of an- 
gles and converted to q-space for sample 2 in Figure [3] 
(black lines), and the azimuthal average of the SAXS 
pattern for a similar sample (red dots). The spectra 
have been vertically offset for clarity. Inset: Posi- 
tion of the peak as a function of angle. Schematic: 
the source and sample are fixed and the detector is 
rotated. 



kapton tape. The peaks from the optical and SAXS 
measurements match when n e = 1.29, implying that 
<j) s = 0.54±0.05. The higher value of (j) s for the thick 
films relative to the thin films is most likely the re- 
sult of the different quench-rates used to make the 
samples. In the spin casting of thin films, the wa- 
ter evaporates within seconds, while water evaporates 
from the thick films over a few hours, allowing par- 
ticle rearrangement which leads to a higher <j) s . In 
the thick film, we find that £ = 940 nm, only a few 
particle diameters. 

These isotropic films mimic the essential optical 
properties of bird feathers that have structural color 
from isotropic nanostructures. A photograph of an 
example of these feathers from the crown of Lepi- 
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dothrix coronata is shown in Figure [5^i, and a trans- 
mission electron micrograph (TEM) of the color- 
producing structure is shown in Figure^. Normal- 
ized I(q) from both L. coronata and sample 2 from 
Figure [3] are plotted as a function of q/qo in Fig- 
ure]^. The SAXS patterns reveal similar structures 
out to the third peak in I(q). Beyond that, the thick 
film has additional peaks that arise from the unifor- 
mity of the spheres used in the sample. We compare 
the performance of the feathers and films by plot- 
ting optical scattering spectra for both at 20° in Fig- 
ure [5]i. The scattered optical intensity peak for L. 
coronata is narrower than the film: the full-width at 
J(A)//(A ) = 0.90 is 49 nm for the feather and 82 
nm for the film. Similarly, L. coronata has a nar- 
rower first peak in I(q). L. coronata also displays 
less scattering at shorter wavelengths. This may be 
due to a significant difference in the two structures: 
the feathers have spheres of air in a high-index of re- 
fraction background whereas the films have spheres 
of a high-index in a background of air. 

Non-iridescent structural color is realized when an 
isotropic structure has a peak in I(q) and wavelength- 
independent scattering is suppressed. Using a bidis- 
perse mixture of spheres produces isotropic struc- 
tures. Controlling film thickness or absorption length 
are effective at enhancing contrast in reflectance spec- 
tra. These films may have applications in coatings, 
cosmetics and textiles. The basic optical properties of 
structurally-colored feathers have been reproduced, 
but more work must be done before the biomimetic 
samples perform as well as their biological counter- 
parts. 
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Figure 5: Comparison with feathers, a) Photograph 
of crown feathers from L. coronata, the field of view 
is 1 cm. b) TEM image of air spheres in beta- keratin 
from L. coronata, the field of view is 5.5 jam. c) 
Azimuthal averages from SAXS measurements of L. 
coronata (dashed line) and a thick film with the same 
[CB] as sample 2 in Figure [3] (solid line). Inset: SAXS 
pattern from L. coronata, the field of view is 0.14 
nm -1 and the gray values are logarithmic in intensity. 
d) Comparison of optical scattering data taken at an 
angle of 20° (experiment geometry shown in Figure |4| 
from L. coronata and sample 2 from Figure [3] The 
line colors correspond to the apparent color of the 
feathers and the thick film. 



5 



1 Experimental 

Monodisperse PS spheres were synthesized using 
a surfactant-free polymerization technique. [19] 
Particles were electrostatically stabilized by co- 
polymerization with sodium 4-vinylbenzenesulfonate. 
After synthesis, particle suspensions were washed by 
centrifugation and resuspension at least three times 
with DI water. After washing, all particle suspen- 
sions were adjusted to (j)ps = 0.3. Particle sizes were 
determined by SEM image analysis using a Philips 
XL-30 ESEM with an accelerating voltage of 10 kV, 
after being coated with a thin layer of gold. 

To prepare bidisperse suspensions, equal volumes 
of two monodisperse suspensions were mixed by 
pipetting approximately 20 times and vortexing for 
at least 30 seconds. Prior to spin casting for thin 
films or water evaporation for thick films, the sus- 
pensions were sonicated for at least 20 minutes. Thin 
films were spin cast onto glass coverslips cleaned with 
ethanol using a Headway Research, Inc. PWM32-PS- 
R790 spin coater. Typical spin speeds were between 
500 and 5000 RPM, the spin speed determined the 
final thickness of the film. [20] 

Cabot Vulcan XC72R GP-3919 carbon black was 
suspended in DI water at [CB] = 4.2 wt% with 1 wt% 
Pluronic F108 to stabilize the CB. Different volumes 
of the CB suspension were added to bidisperse sus- 
pensions of spheres to produce samples with different 
final [CB]. 

To estimate the extinction lengths quoted in Ta- 
ble [T] the transmission spectra of aqueous suspen- 
sions of CB were measured using an Ocean Op- 
tics USB 650 Red Tide spectrometer. Suspensions 
ranging from 8xl0 -5 to 2xl0 -2 wt% were used. 
Measurements were made with three different path 
lengths: 1 cm, 0.04 cm, and 0.03 cm. The ex- 
tinction length for all path lengths was proportional 
to [CB] -1 , allowing us to extrapolate the extinction 
length to the [CB] range used in our thick film sam- 
ples. We use extinction instead of absorption because 
our measurement did not discriminate between ab- 
sorbed and scattered light. 

SAXS measurements of biomimetic samples were 
carried out at beamline X9 at the NSLS, Brookhaven 
Naitonal Laboratory, using a Rayonix Mar 165 CCD 



detector. The X-ray energy was 7 keV (a wavelength 
of 1.771 A) and the sample-to-detector distance was 
5 m. The conversion from the detector image to recip- 
rocal space was calibrated using the diffraction pat- 
tern from a standard silver behenate sample (qo = 
0.1076 A -1 ). The first ring from silver behenate is 
out of the angular range covered by the detector at 
7 keV, therefore the standard pattern was collected 
with the same scattering geometry but at an X-ray 
energv of 15.65 keV (0.792 A). The finite beam size 
at the detector corresponds to FWHM q resolution 
of 0.00045 A -1 , which is -17% of the width of the 
fringe produced by the spheres (7r/120 nm « 0.0026 
A -1 ). Biomimetic samples for SAXS measurements 
were prepared on 0.0025" -thick Kapton Tape, pur- 
chased from McMaster-Carr (catalog no. 7648 A33). 

SAXS measurements of the feathers were carried 
out at beamline 8-ID-I at the APS, Argonne National 
Labs, as described in Dufresne et al. [4 

To simulate the structure of our samples, we cre- 
ated mechanically stable packings of bidisperse fric- 
tionless spheres in cubic cells with periodic bound- 
ary conditions using a simulation protocol in which 
we successively compress/decompress soft particles 
and then apply conjugate gradient energy minimiza- 
tion until particles are just at contact (up to a pre- 
scribed energy threshold). [21] The energy thresh- 
old, initial packing fraction (/>o, and increment in <j) 
were set to 10 -8 , 0.2 and 10 -4 , respectively. We 
selected particle size and number ratios to match 
the experiments. From 50 independent initial con- 
figurations, we obtained an average packing fraction 
= 0.63, which is relatively insensitive to the spe- 
cific parameters of the packing-generation protocol. 
From the particle centers, we calculated the partial, 
Su(q), S ss (q), and Si s (q), and total structure factor 
S(q) = xiSu(q) +x 3 S 33 (q) + 2y/xix 3 Si 3 (q), where l,s 
signify large or small particles, x is the number frac- 
tion of the indicated particle, and q is the wavevector. 
The partial structure factors and the theoretical form 
factors for each sphere size were used to calculate 
the total scattering intensity I(q) oc xiSu{q)Ff{q) + 
2^xWsSi s (q)F l (q)F s (q) + x s S ss (q)F 2 s (q). 

Optical transmission spectra were performed using 
a Hitachi U-2001 spectrophotometer. Optical reflec- 
tion and scattering spectra were performed with a 
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custom-built setup in Hui Cao's laboratory described 
in detail in Noh et al. [T8] 
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